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Abstract 

The conservative dynamics of two point masses given in harmonic coordinates up to the third 
post-Newtonian (3pN) order is treated within the framework of constrained canonical dynamics. A 
representation of the approximate Poincare algebra is constructed with the aid of Dirac brackets. 
Uniqueness of the generators of the Poincare group resp. the integrals of motion is achieved by 
imposing their action on the point mass coordinates to be identical with that of the usual infinites- 
imal Poincare transformations. The second post-Coulombian approximation to the dynamics of 
two point charges as predicted by Feynman- Wheeler electrodynamics in Lorentz gauge is treated 
similarly. 
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I. INTRODUCTION 



High post-Newtonian accurate description of general relativistic dynamics of compact 
binaries in harmonic gauge has many applications in relativistic astrophysics, notably in 
connection with binary pulsars and future gravitational wave astronomy, see e.g. [1-3]. Inspi- 
ralling compact binaries are even the most promising sources to be detected by ground-based 
interferometers such as LIGO, VIRGO and GEO600. The corresponding higher order post- 
Coulombian approximation offers a simpler analogue to post-Newtonian dynamics relevant 
also on its own [4, 5]. 

The approximate analytical dynamics of compact binary systems in general relativity is 
most often treated in harmonic coordinates (see [6] and references therein). Quite recently 
the dynamics of binary point masses has been completed to the third post-Newtonian (3pN) 
order [7-13]. Hereby results derived by means of the canonical formalism of Arnowitt, 
Deser and Misner have been confirmed [14-20]. In approaches based on the use of harmonic 
coordinates, the dynamics was first obtained under the form of ordinary second order 3pN 
equations of motion satisfied by the particle trajectories. The Lagrangian corresponding to 
the conservative part of the motion turns out to be of higher order in the time derivatives 
of the point mass coordinates [21]. This feature is shared by the Lagrangian of Feynman- 
Wheeler electrodynamics in Lorentz gauge derived by Kerner [22, 23]. In both cases, Euler- 
Lagrange equations, of third order or higher, admit a wide class of solutions including 
physically irrelevant ones that do not reduce to the Newtonian solution in the limit where the 
speed of light c tends to infinity. This can be seen directly from the number of independent 
initial data including higher than first order time derivatives of the position variable. When 
demanding the Newtonian limit, it becomes possible to re-derive the ordinary second order 
equations of motion by an iterative order reduction procedure. 

The higher order property of these Lagrangians actually arises from the fact that the 
dynamics in harmonic coordinates as well as in Lorentz gauge for electromagnetism are 
approximately Poincare invariant. Indeed, the so-called no-interaction theorem by Martin 
and Sanz [24] states that Lagrangians of point particles derived by means of a slow motion 
approximation from some classical field theory must contain higher order derivatives from 
second order level in the 1/c 2 expansion, if approximate manifest Poincare invariance is 
maintained. For arbitrary approximately Poincare invariant point-particle dynamics, higher 
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order derivatives must be contained in the Lagrangian only from the third order in powers 
of 1/c 2 [25]. It can also be shown that, if exact Poincare invariance of a system with finitely 
many degrees of freedom is required, (i) Lagrangians including interactions must depend 
on time derivatives of infinite order [26], (ii) for point particle systems, the positions may 
not be chosen as canonical coordinates in Hamiltonian formalism (which reflects the time 
non-locality due to retardation) [27-29]. 

For the approximate dynamics, appropriate contact transformations lead to an ordinary 
Lagrangian but in a non-harmonic grid [30] resp. in a non-Lorentzian gauge. The trans- 
formed representation of the dynamics can be described by means of an ordinary canonical 
formalism. Another approach consists in constructing the canonical formalism corresponding 
to the dynamics directly in the original frame. 

The original higher order Lagrangian is of singular type because of the higher order 
derivatives occurring in "small" corrections i.e. in terms of higher order in powers of 1/c 2 . 
The Hessian is thus multiplied with some positive power of 1/c 2 . Therefore, it is non- 
invertible on R[l/c 2 ]/(l/c 2n+2 ), the ring of real polynomials in 1/c 2 modulo C(l/c 2n+2 ) we 
are working on at the nth order level of approximation. Indeed, in expressing the highest 
order derivative as a function of the others, the Euler-Lagrange equations are multiplied with 
the inverse matrix of the Hessian, so that, in particular, the Newtonian part is multiplied 
with some power of c 2 and the Newtonian limit does not exist anymore. Independently from 
the l/c 2 -power prefactor, what we shall call the "matrix part" of the Hessian may be singular 
by itself, leading to an additional singular structure. 

The canonical formalism for singular Lagrangians goes back to Dirac [31, 32], as well as 
Anderson and Bergmann [33]. The first canonical treatment of dynamics derived from a slow 
motion approximation of a classical relativistic theory described by a singular Lagrangian of 
higher order (in the time derivatives) is due to Jaen, Llosa and Molina [34]. Specializing on 
a class of approximate Lagrangians of a certain structure in the l/c 2 -power expansion and 
demanding an invertible "matrix part" of the Hessian, they developed a method aimed at de- 
riving an explicit expression for the Hamiltonian as well as the Dirac brackets. They applied 
their formalism to the 2pC dynamics in Lorentz gauge. However, the resulting Hamiltonian 
was not correct because of computational errors. The first correct 2pC Hamiltonian is due to 
Damour and Schafer [35], their approach having been detailed in [30]. Later, Saito, Sugano, 
Ohta and Kimura proposed a method how to treat general higher order singular Lagrangians 
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in canonical formalism. They proved the equivalence of Lagrangian and Hamiltonian formu- 
lations for singular Lagrangians of higher order [36, 37]. A similar analysis has later been 
performed by Gracia, Pons and Roman-Roy in a geometrical framework [38]. The formal- 
ism given in [36] was applied to a class of 2pN Lagrangians to which the post-Newtonian 
Lagrangians in harmonic coordinates do not belong. It was also used by Ohta and Kimura 
for investigating aspects of 2pC Feynman- Wheeler dynamics in Lorentz gauge [39]. Note 
that the approaches of articles [36] and [34] are crucially different. The singularity arising 
from the fact that we are working on the ring M[l/c 2 ]/(l/c 2n+2 ) is indeed not considered in 
Ref. [36]. 

The aim of this paper is to formulate the canonical formalism for the conservative part of 
3pN dynamics in harmonic coordinates as well as for 2pC Feynman- Wheeler electrodynamics 
for two particles in Lorentz gauge and to analyze the dynamics in this framework. For the 
formulation, we use a similar method as the one developed in Ref. [34] generalized to 3pN 
conservative binary dynamics. For the first time we give the 3pN Hamiltonian in harmonic 
coordinates, the corresponding Dirac brackets, as well as a canonical representation of the 
Poincare algebra of the 3pN and 2pC dynamics. 

The canonical description is always helpful for a better understanding of the dynamics. 
It is an extremely elegant tool to derive features such as symmetries and integrals of motion. 
The latter quantities, computed in harmonic coordinates at the 3pN order and specialized to 
the center of mass frame [40], are useful for the description of inspiralling compact binaries 
relevant as sources of gravitational waves. They allow the derivation of an analytic para- 
metric "generalized quasi-Keplerian" solution to the 3pN accurate conservative equations of 
motion for compact binaries moving in eccentric orbits [41]. This is relevant in particular 
to construct post-Newtonian search templates for the detection of gravitational waves or to 
compare the numerical and post-Newtonian descriptions of such systems. Our integrals of 
motion prove to be consistent with those computed by Andrade, Blanchet and Faye [21] in 
Lagrangian formalism, providing a powerful cross check for the results. 

Furthermore the inclusion of spin in post-Newtonian binary dynamics using covariant 
spin supplementary condition also results in a dynamics that is described by a higher order 
singular Lagrangian or Hamiltonian when staying in harmonic coordinates [42]. Especially, 
the investigation of its canonical description derived using the methods in the present article 
will likely be useful for the prediction of gravitational wave templates. 
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This article continues work initiated by Stachel and Havas who derived in 1976 the 
Hamiltonians describing a class of dynamics including the lpN and lpC ones and who 
computed the integrals of motion corresponding to the approximate Poincare invariance 
[43]. They announced a further article, where special interactions allowing the choice of the 
spatial coordinates as canonical coordinates were to be treated up to second order, but this 
article was never published. We do not follow the program they had initially designed but 
rather concentrate on physically relevant interactions incompatible with the latter choice of 
canonical coordinates. 

The plan of the paper is as follows. In Sec. II, we outline the general constrained La- 
grangian formalism for Lagrangians containing higher order time derivatives. We also show 
how to derive a full time-stable set of Lagrangian constraints not only for Lagrangians having 
a similar structure as in [34], but also for cases where the "matrix part" of the Hessian is not 
invertible. In Sec. Ill, we outline the theory of the corresponding Hamiltonian formalism. 
Sec. IV is dedicated to a short description of the Poincare algebra and its action on spatial 
coordinates. In Sec. V we apply the preceding results to the 3pN dynamics of two point 
masses, showing the explicit Poincare invariance and deriving the corresponding integrals of 
motion. The 2pC Lagrangian of Feynman- Wheeler electrodynamics for two charged point 
masses is treated similarly in Sec. VI. Finally, in Sec. VII, we summarize and discuss our 
results. 



II. HIGHER ORDER SINGULAR LAGRANGIAN POINT MASS DYNAMICS 

We start from the action integral of a higher order Lagrangian L that does not depend 
explicitly on time. It simply reads [44, 45] 

S = dtL(q,qW,...,qW), (1) 
J t 

where q is a short notation for the set of / independent variables {q^}, \i = 1, ...,/, and 
where denotes the set {q$} of their ith derivatives with respect to time t. The highest 
order of time derivative appearing in L is denoted by n. From the action principle 8S = 0, 
we draw the generalized Euler-Lagrange equations, 
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with (j, — 1, . . . , / , d/dt denoting the total time derivative. By collecting the terms of the 
/xth equation that do not depend on into a single function and isolating the highest 
order derivatives q( 2n \ we may rewrite Eq. (2) as 
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This relation shows that the highest order time derivatives always occur multiplied with the 
Hessian matrix 



+ /,(g,...,g( 2 "- 1 )) = 0. 
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so that the Euler-Lagrange equations can be solved for q^ as a function of the configuration 

space variables q^ 2n ~^ if and only if the Hessian is invertible. 

We now want to specialize to Lagrangians derived within slow-motion approximation 
schemes of relativistic theories. These are for instance the Lagrangian describing the con- 
servative part of post-Newtonian dynamics that has been determined up to third post- 
Newtonian order in Ref. [21], or the Lagrangian describing the nth post-Coulombian dy- 
namics of two particles in the Feynman- Wheeler theory of electromagnetism [23]: 
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where the letter x is used to refer to configuration space variables in order to emphasize 
the physical interpretation of x^ as the jth derivative of the ith component (i = 1,2,3) 
of the position vector of the particle a = 1,2; where Xo' stands for the jth time derivative 
of the position vector of particle a, and x^ = {ia '} for the set of all configuration space 
variables that are time derivatives of jth order (x = x^). Round brackets (. .) indicate 
the scalar product, e.g. (x^x^) = Yli=i x ii x 2i '■> if both vectors are identical we denote, 
e.g. (x^) 2 . We further introduced the operator D a , which represents the time derivative 
acting exclusively on the configuration space variables belonging to the particle numbered 
a. r is the absolute value of the relative separation vector, and e a denotes the charge of the 
particle a. Note that, since the square root may be expanded into a binomial series up to 
order l/c 2n+2 , it makes sense to consider iV-particle Lagrangians of the type [34] 

1 N n 



a=l 



s=0 



with e = js. Whereas for the 3pN dynamics there is actually no a^ 3 ) dependence, the 
formalism can be adapted to this case (see Sec. V). The Hessian of the above Lagrangian 

TT _ r n 9 2 V n 

n-aibj - t , s , , {I) 
dx ai dx bj 

is not invertible on the ring R [e] / (e n+1 ) on which the approximation scheme is defined. 
Thus, the Lagrangian is singular and the system is subject to constraints which are now to 
be determined. 

The primary Lagrangian constraints are given by all independent linear combinations 
of the Euler-Lagrange equations that do not contain 2nth order time derivatives, imposing 
thereby conditions on the configuration space variables (cf. e.g. [46, 47]). From Eq. (3) we 
see that they can be derived by contracting the Euler-Lagrange equations with some null 
vectors of the Hessian (7), we shall refer to as X r . Let us first suppose that the "matrix part" 
— M Vn (n ) of H a ibj is invertible [341, a restriction we will skip later. Then, in our approximation 

dx ai dx bj 

scheme, the null vectors of the Hessian are those that are multiples of e. The contraction of 
the Euler-Lagrange equations with the canonical basis vectors of R 37V multiplied by e yields 
all the primary constraints. In the notation of Eq. (3) with the generalized coordinates being 
the spatial coordinates x and their time derivatives, these are given by 

ef bj (x,...,x( 2n -V) = 0(e n+1 ). (8) 

Since — ( ^ y " (n) is invertible, there are no more independent ones. Requiring the special form 
(6) of the Lagrangian, the Euler-Lagrange equations read 

n 

tX W + J2e s A as (x, ..,x (2 ^) = 0{e n+1 ), (9) 



-»>r. 

s=0 



dY dV a 



r=0 



dtj dx^ 



(r)> 



with a = ai and m a = m a . Now, by means of Eq. (9), we rewrite the primary constraints 
(8) as 

n—l 

em a x^ =eY,e s A as (x,...,x^) + 0(e n+1 ). (10) 

s=0 

We observe that they coincide up to the factor e with the equations of motion of lower order 
in e. Starting with the above equation, we can derive a minimal stable set of constraints as 
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explained in appendix A, 

1 



T (2+r) _ _ 



+ 0(e n+1 ), (11) 



J2 £SB *,2+r,s( X , x{1) ) 
s=0 

for r = 0, ...,2n — 3. The precise definition of the functions B a2 + r ,s from the A QS 's is 
specified in the appendix. It is worth noticing that the constraints corresponding to the case 
r = agree with the equations of motion after they have been iteratively reduced to order 
two in the time derivatives by removing higher order time derivatives with the help of the 
equations of motion of lower order in e. Similarly, the additional constraints agree with the 
appropriately reduced time derivatives of the reduced equations of motion. This justifies the 
preceding statement saying that the constraints emerge by requiring the Newtonian limit. 



The "matrix part" — ( n ) Vn ( n ) of the Hessian of post-Newtonian Lagrangians linear in the 

ax bj 

accelerations, is not invertible. This can be cured by adding so-called double-zeros. While 
this may change the rank of the "matrix part" of the Hessian, it does not influence the 
order-reduced equations of motion [48, 49]. Because of the agreement between the latter 
equations and the constraints (11), established for an invertible matrix — ( ^ v " (ra) , we do not 
expect double-zeros to change the constraints either. A closer investigation shows that this 
is indeed the case. We can state even more generally that the expression for the constraints 
and the construction of the Hamiltonian remains unchanged if the "matrix-part" of the 
Hessian is non-invertible. 



We thus suppose that 



d 2 Vu 
dx^dx^ 

at bj 



is not invertible and has a constant rank 3iV — R all 



over the configuration space. Then, in the considered approximation scheme, all multiples 
of £ are still null vectors of the Hessian, but there also R additional null vectors say Ap Q , 
p — 1, R, of order zero in e. The primary constraints are obtained by contracting any of 
them with the equations of motion. We shall first consider the constraints emerging from the 
contraction of the Euler-Lagrange equations with the canonical basis vectors of M 37V times 



e. Since the regularity of 



d 2 V n 



dx<$dx§ 



- is actually not used in the derivation of the minimal 



time stable set they belong to, this set is again given by Eq. (11). The primary constraints 
generated by the additional null vectors \ pa read 



s=0 



=0(e n+1 ), 



(12) 
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where p = 1,...,R and where J2 a denotes the sum over all pairs a = ai with a = 1,2 and 
i = 1,2,3; \ pa may depend on x,...,x^ for the Hessian itself possibly depends on these 
variables. We must now examine the additional restrictions imposed by Eqs. (12) to the 
already derived constraint surface given by the set of relations (11). On this surface, by 
definition all higher order derivatives x^, ...,x^ 2n ~^ entering Eqs. (12) may be expressed by 
means of the coordinates x, x^ with the help of Eqs. (11). We find thus 



^ ^ Ap a (a;, . .., x^ ^ 

a 

Evo^ (1) ) 



a%a ^ ~t~ ^ ^ ^ A as (x, . . . , X^ ^ ) 



-m 

n 



+ 0(e n+1 ) 



+ 0(e n+1 ) « 0(e n+1 ). (13) 
(ii) 



-m a x^ + S £ j e s B aAs (x,x {l) ) 

We emphasize the fact that these relations only hold on the constraint surface by using the 
weak equality symbol The system of equations (13) tells us that the constraints resulting 
from the additional null vectors of the Hessian are already fulfilled on the surface defined by 
the constraints (11). It is satisfied for all times due to the time-stability property, so that 
the seemingly additional constraints are covered by the set (11). In short, the additional 
null vectors do not generate additional constraints. This fact enables us to perform the 
transition to Hamiltonian formalism regardless of the invertibility of the "matrix part" of 
the Hessian. Moreover, the Lagrangian constraints of the considered dynamics are still 
identical with the reduced equations of motion or their reduced derivatives. In particular, 
double zero terms, though they may change the rank of the Hessian, do not influence the 
formalism as long as the general structure (6) is maintained. We observe that according 
to above computation, unlike the usual theory, there are no arbitrary functions of time 
emerging in the dynamics, even if the contraction of some of the additional null vectors with 
the Euler-Lagrange equations vanishes identically [46, 47]. This is ultimately a consequence 
of the linear independence of the equations of motion at lowest order. 

III. HIGHER ORDER SINGULAR CANONICAL FORMALISM 

A system described by a Lagrangian of higher order allows for a canonical description 
with phase space variables and canonically conjugate momenta with j = 0, ...,n — 1 
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[50, 51]; II, is the set of jth so-called Ostrogradski momenta 11,1, ■ ■■,'Kjf defined by 

n » ^ £ N (14) 

with j = 0,...,n— l,/i = 1, ...,/, n being the highest order derivative entering the Lagrangian 
and / the number of degrees of freedom. For j — 0, n — 2 we may write alternatively 

/ 

n*. = E(- 1 ) n " j_1 ^ 2n " j_1) ^ + ? (an - J "- 2) ) (is) 

showing that the highest order time derivatives occur multiplied with the Hessian. For 
j = n — 1 we have 

dL . . 

It is of the form (15) with Kj^ = if L is quadratic in the qjfi. (This special case has 
not been accounted for in Ref. [34].) Let us first assume that the Hessian is regular. In 
this case, the Ostrogradski transformation can be inverted by using an iterative algorithm. 
The implicit equation (16) is locally solvable for the and yields q^ n \q, <2^ n-1 \ n n _i). 
Having computed the variables g( n+l ) with i < j, we can invert the equation for the n— j — 1th 
Ostrogradski momentum for q^ n+ ^(q, q^ n ~ l \ n n _j_i, II n _i). 
The Hamiltonian of the system is 

H = - L(q, .., q^\ q^(q, g^ 1 ), H^)) (17) 

n-1 f f 

+ EE^? +1) +E n Mrf ) (?' ....(z^.n^o, 

while the Hamiltonian equations of motion take the familiar form 

HP = «v (18a) 
s n - ^ (18b » 

with fc = 0, ...,n — 1. Introducing the (Ostrogradski-) Poisson bracket 

w^J^T——-—— (i9) 

the time evolution equations for a smooth function F of the phase-space variables and time 
t takes the form 

d f)F 

- F={FlBi + -. (20) 
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Let us now turn to the case where the Hessian is singular. Then, the Ostrogradski trans- 
formation (q, g( 2n_1 )) — > (g, n , II n _i) is not invertible anymore, or equiva- 
lently, the phase space variables considered as functions of q, ...,g( 2n_1 ) are not all indepen- 
dent. If the Hessian has rank f—r (for a maximum rank /), they are linked by r independent 
relations. These are the primary constraints 

^ a (q,...,q^ n - 1 \U ,...,U n . l )=0 (21) 

following from the definition (14) of the momenta. To make sure that the resulting constraint 
surface T be a submanifold of phase space, we impose the "regularity condition", demanding 
that zero be a regular value of ^ regarded as a map on phase space to W. 
The Hamiltonian if as a function of the configuration space variables reads 

n-l / 

H(q, q^)) = - L(q, q^) + £ £ H^, q^'^q^. (22) 

j=0 m=0 

The remarkable fact is that the coordinates q^ n \ q( 2n ~V appear only through the com- 
binations Iijn(q, g( 2ri-J ' -1 )) due to the particular form of the dependence of L and n jM 
on those coordinates. Hence, H actually depends only on q,...,q( n ~^ and n , ...,n n _i. 
This can be verified in a similar way as in the absence of higher order derivatives [52, 53]. 
Therefore, we may view the Hamiltonian as a function of the phase space variables 
H = H(q, q^ n ~ l \ n , ...,n n _!) although it is not unique in the case where the Hessian is 
not invertible but defined modulo a linear combination of the primary constraints YT a =i ° a ^a, 
with c a being functions of the phase-space coordinates [52, 53]. The time evolution of a 
smooth function F of the phase space coordinates and time is given by 

j t F = {F,H} + J2u a {F,* a } + ^, (23) 

a=l 

where u a , a — 1, ...,r are extra parameters, and { . , . } refers to the Poisson bracket (19). 
Additionally, at this level, ^ a = has to be imposed on the motion. The relation d^f a /dt = 0, 
a = l,...,r may entail restrictions to the u a, s and/or lead to new, secondary, constraints. 
Using again the time stability property, we may get further conditions, and so on. At the 
end we are left with a complete set of, say, K, time stable constraints including the primary 
ones 

tt fc (g, n , n n _o = o, k = 1, .., k, (24) 
11 



subject to the same regularity assumptions. 

In constrained dynamics, there is an important distinction between two types of con- 
straints: (i) First class constraints are characterized by the property that their Poisson 
brackets (19) with all the other constraints vanish on the constraint surface, (ii) second-class 
constraints have at least one non- vanishing Poisson bracket on T. If the Poisson-bracket ma- 
trix D, defined in components as D^i = {^k, ^1} for k, I — 1, K, has rank A on T, then the 
set of constraints (24) can always be linearly transformed into an equivalent set consisting of 
K — A first-class constraints and A second-class constraints with invertible Poisson-bracket 
matrix on Y [52]. Conversely, if no first class constraints can be obtained from Eqs. (24) by 
linear transformations, the Poisson bracket-matrix D is invertible on T. In this case, the 
Dirac bracket of two phase space functions may be defined by 

K 

{F, G}* = {F, G}-J2 {F, ^k}D^{%, G}, (25) 
k,i=i 

where D~ l is the inverse matrix of D, existing a least in a neighborhood of T. 

The Dirac bracket keeps important features of the Poisson bracket; namely it is bilinear, 
antisymmetric, acts as a derivation on each argument, and fulfills the Jacobi identity. It 
actually is the restriction of the Poisson bracket on the constraint surface [54]. Finally, it 
has by construction two additional important properties 

{# fe ,F}* = 0, (26) 
{G, F}* m {G, F}, (27) 

valid for arbitrary phase-space functions F and for functions G of first class, i.e. functions 
whose Poisson bracket with any constraint vanishes on T; "rs" represents the weak equality 
holding only on T. We are thus allowed to simplify expressions entering the Dirac bracket 
by using the constraint equations, which amounts to setting the \I/fc to zero, before the final 
computation of the bracket. With the help of such a tool we are in position to reformulate 
the time evolution equations equivalent to Eq. (23) on the constraint surface as 

_ F=Wfll} . + _, ( 28 ) 

where H r is the reduced Hamiltonian, derived from H = —L + Y^jZo S^=o ^>9m by 
eliminating all the other coordinates in favor of the coordinates of the constraint surface 
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with the help of the system's constraints [34, 47, 52, 53]. The preceding canonical formalism 
has been shown to be equivalent to the corresponding Lagrangian formalism in Ref. [36, 38]. 

We now want to specialize on systems described by a Lagrangian of the form (6). From 
Eq. (14) we compute the Ostrogradski momenta, 

n ja = m a x^5 j0 + e j+1 $ ja (x, x^'^) + 0(e n+1 ), 

s=0 1=0 v 7 ux a 

The Lagrangian constraints (11) are already known. Owing to the equivalence theorems 
proved in [36, 38], they can be translated into constraints of the canonical formalism by 
keeping the first n — 2 relations and eliminating x^ 2 \ x( 2n_1 ) from the n identities (29) by 
means of Eqs. (11) (cf. Ref. [34], in which this method is applied whereas the equivalence 
had not yet been formally stated). We find 

1 " 

u ra = a£) £SB °,r,s(x, = 0(e n+1 ) (30a) 

1710 3 = 

uj la = x& - — [n 0Q - e$ 0a (x, x^)] = 0(e n+1 ) (30b) 

TYl a 

Xja = n ja -ei +1 $ ja (x,xW) = 0(e n+1 ), (30c) 



with r = 2,.. .,n — 1 and j = l,...,n — 1. The &j a (x, iW)'s are derived from the 
$j a (x, x( 2n ~i~ r> } 's by eliminating higher order derivatives with the help of Lagrangian 
constraints. As the constraints (30a-30c) obviously fulfill the regularity conditions, they 
define a constraint surface Y that is a submanifold of phase space. As coordinates, we may 
choose e.g. x and n , another possibility would be x and x^ l \ 

First class Hamiltonian constraints occur only if there are configuration space variables 
that are arbitrary functions of time in Lagrangian formalism. We have already seen that 
these are not present in our case so that there are no first class constraints in the dynamics 
under consideration. This result can also be proven by direct computation of the Poisson- 
bracket matrix (cf. appendix C and [34]). It can be checked that 

^ _ l{Xka,Xrf3} {Xka,Urf3}\ 
\{Uka, Xr@} {^fco, Urp} I 

is indeed invertible, which implies in particular that the number of constraints is even. Hence 
the constraint surface T has an even dimension, say g. Now, by an iterative procedure, D^ 1 
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may be expressed by means of submatrices of D (cf. appendix C and [34]). Adopting the 
splitting 

D " = (i ■ 1321 

we may then write the Dirac bracket of two functions / and g on the phase space as 

n—1 n—1 

{f,g}* = {f,g}~ Y ^2{fiXka}X ka:lf3 {xif3,g} - Y ^2{fiXka}Y kat ip{uip,g} 

k,l=l a,0 k,l=l a,P 

n—1 n—1 

+ Y ^2{f^ U ha}yip,ka{XW,9} ~ Y 53{/> W *a}^k*,J/3{wj/9,0}- (33) 

k,l=l a,f3 k,l=l a,f) 

This bracket defines a symplectic form on the constraint surface. According to the theorem 
of Darboux, we can locally find canonically conjugate coordinates Q, P of T such that it 
takes the familiar shape 

{F ' = S ~ dQ a dP a (34) 

We shall see in Sec. V and Sec. VI that, for the dynamics investigated in this article, it 
will even be possible to exhibit global coordinates of this kind, which will greatly simplify 
computations involving Dirac brackets. The global existence of canonically conjugate co- 
ordinates entails that the group of generalized canonical transformations, i.e. the group of 
transformations that leave the Dirac bracket invariant, is necessarily a subgroup of the group 
of canonical transformations on the phase space [47]. 



IV. POINCARE ALGEBRA 

The symplectic structure on the phase space (resp. on the constraint surface) allows us to 
endow the vector space of scalar fields defined on T with the structure of a Lie algebra. The 
Poisson-bracket (resp. Dirac-bracket) relations satisfied by the generators of the infinitesimal 
(generalized) canonical transformations corresponding to the action of a transformation 
group are known to be identical to the Lie-bracket relations of the generators of this group. 
In other words, there exists a Lie algebra homomorphism between the Lie algebra of the 
transformation group and that of the generators of the infinitesimal (generalized) canonical 
transformations (see e.g. [55]). 
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Knowing the latter Lie algebra, we are enabled to reconstruct the symmetry group locally. 
This can be done globally only for simply connected groups (unlike the Poincare group). 
As we aim at determining the symmetry group of some system in the phase space or on 
the constraint surface locally, it will be sufficient for us to consider the Lie algebra of the 
generators of infinitesimal (generalized) canonical transformations (see e.g. [56]). 

By construction, the conservative pN dynamics in harmonic coordinates as well as the 
pC dynamics of the Feynman- Wheeler theory in Lorentz gauge are approximately Poincare 
invariant. Therefore, in a neighborhood of the identity there exists an approximate repre- 
sentation of the Poincare group as a generalized canonical transformation group on T. This 
means that on the constraint surface, there are generators H, P iy Ji, d with i — 1,2,3 of 
infinitesimal generalized canonical transformations that approximately fulfill the Poincare 
algebra with respect to Dirac bracket: We have thus 

{Pi'VY = o(^), (35a) 

3 



c 2(n+l) 

ft, J jY = E ^ + ° ( ^Ld ) ' ( 35b ) 



k=i 

3 



{J,, Jj}* = Sikh + O (— ^ ) , (35c) 
fe=i ^ ' 

{ H > p iY = (^TT))' ( 35d ) 

{^ J iY = ( 35e ) 



c 2(n+l) 
3 



1 / 1 

{Gi,Gj}* = --J2 e vkJk + I 

fc=l ^ 



(35f) 



{ Gl ,HY = P i + o(^^y (35g) 

{J t , Gj}* = J2 SikGk + O (-^) , (35h) 
fc=i ^ ' 

{G t ,P 3 Y = ^H + o(^^j, (35i) 

where i,j = 1, 2, 3 are spatial indices and where n is the order of approximation (n — 3 for 
the 3pN, n = 2 for the 2pC dynamics). In the non-relativistic limit c — > oo, the algebra 
(35a-35i) reduces to the Galilean one. In particular Eq. (35i) becomes 

{Gi, PjY — &ijM, (36) 
15 



M being the total non-relativistic mass of the system. Since all quantities are defined on 
the ring M[l/c 2 ]/ (l/c 2(n+1 )), strictly speaking, the term Mc 2 appearing in the generator 
H is not allowed. However, because it is a mere constant, it does not change the dynamics 
of the system. We may keep it for it suggests the physical interpretation of H as the total 
conserved energy. To harmonize the argument, we shall work with a slightly modified form 
of Eq. (35i), namely 

{ Gl , P 3 Y = 5 13 {M + ±H) + (-j-Lj) . (35i') 

The statement that a dynamics is Poincare invariant usually signifies more than the 
existence of generators satisfying the Poincare algebra. Namely, that the generators of the 
infinitesimal canonical invariance transformations H, Pi, Ji, Gi are the generators of the linear 
representation of infinitesimal time translations, spatial translations, spatial rotations and 
Lorentz boosts acting on the particle components. In short, the members of the Poincare 
group are to have the usual action on the physical system, when the phase space coordinates 
x(t) are interpreted as the positions of the particles in Cartesian coordinates parameterized 
by the time t. This is often emphasized by referring to this type of dynamics as "manifestly" 
Poincare invariant. 

In this perspective, we shall investigate one by one the actions of the preceding gener- 
ators on the phase space variables, acting on the coordinates x(t) regarded as the spatial 
coordinates of the system of N point particles at time t. The time evolution of x^(t), 11^ (t) , 
j = 1, ...,n — 1 is given by the trajectories on the constraint surface; the argument t will be 
dropped below for the sake of simplicity. Throughout this article we shall adopt the point 
of view of active transformations. 

We start with the generator of infinitesimal time translations H. By definition it generates 
a transformation that can be interpreted, when regarded as active, as a translation of the 
particles along their trajectory. Now, the effect of an infinitesimal time translation about r 
is 

r^S = 6rx$ = T{x2\H}- + O , (37a) 

rj t U kai = 5 T U kai = r{U kai , H}* + O (^i) ) • (37b) 

Note the use of Dirac brackets due to the fact that the trajectories are restricted to T. 
From the equations of motion (28), we see that H may be chosen as the Hamiltonian or 
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the reduced Hamiltonian of the system as has been anticipated by the notation. Since H 
can be interpreted as the total energy of the system, the identity {H, H}* = traduces the 
conservation of energy in time. 

The generators Pj and J«, interpreted as linear and angular momentum, are the genera- 
tors of infinitesimal spatial translations and rotations. Even in absence of information on the 
dynamics, we know their Poisson bracket with arbitrary scalar fields defined on the phase 
space. We know furthermore that generalized velocities are contravariant while generalized 
momenta are covariant vectors. Therefore, under an infinitesimal translation about gen- 
erated by Pi, and an infinitesimal rotation about ipi, generated by J«, i = 1,2,3, the phase 
space coordinates representing the particle components transform as 

3 

3=1 
3=1 

3 

o 5 t n m a i ^ ^ € j { n m a i , Pj } , 

3=1 

3 3 

j,k=l 3=1 
3 3 

^ ^ ^ijk^Pj^m ak ^tp^-mai ^ ^ j { ai i 'Jj}j 
j,k=l j=l 

with / = 1, n — 1 and m = 0, n — 1. They lead to the differential equations 

dPi dPi 

8ik8om, — 7— T = 0, (38a) 



o n a (m) 

<9J,; ( m ) OJi 



an ~^2 €i 3 kX aj > ( m ) - ^O^mak, (38b) 



(m) 

j=l ^ajf fc=l 

which fix the momenta P, and Jj up to a constant. The Poisson bracket relations of the 
these generators are well known; we have for instance 

3 



{Pi,Jj} = ^2e ijk P k , (39a) 
k=i 

3 

{Ji,Jj} = J~]e ijk J k . (39b) 



k=i 
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From Eqs. (38a-39b) we derive the unique result 

N 

p, t = ^n 0(M , (40) 

a=l 

n-1 N 3 

J i = ^2Y1 Yl e HkX { ™ ] Tl mak . (41) 

m=0 a=l j,k=l 

If a given dynamics is invariant under spatial translations and rotations, so must be the 
constraint equations and the constraint surface. In this case, the generators of spatial trans- 
lations and rotations are thus first class functions. From Eq. (27), valid on the constraint 
surface, we see then that all Poisson-bracket relations including Pi and Ji also hold on T 
as Dirac-bracket relations. This means that the momenta Pi and Ji displayed above are 
the generators of infinitesimal spatial translations and rotations represented as generalized 
canonical transformations on T. 

The physical interpretation of Gi as generator of infinitesimal Lorentz boosts allows us 
to determine its action on the particle coordinates x ai [43]. An infinitesimal boost about Vj, 
j = 1, 2, 3 acts on the space-time coordinates of a particle following 

1 3 

5t a = — -^VjXgj, (42a) 

° 3=1 

5x ai = -Vit. (42b) 

A particle, located at x ai at time t, is located after the active transformation at position 
x' ai = x a i + 8x a i at time t' a — t + 5t a . In the three-dimensional space, this results in 

x'aiit'a) = X ai{t) + ^X ai = X ai (t) - Vit. (43) 

Let us stress that, because we interpret the transformation as active, we keep the same space- 
time coordinate system and simply boost the particles. We choose the time coordinate t 
to be eliminated in favor of t' a up to first order in Vi. (However, since we are dealing with 
functional identities valid for all the times, we could also proceed by substituting for t.) 
Applying t = t' a + 4, Y? i= 

zl VjX a j(t) and expanding Eq. (43) up to the linear order in Vi, we 

arrive at the relation 

1 3 

<M = *«(0 + ~ 2 J^VjXajiOXaiiO - + 0(v 2 ). (44) 

3=1 
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We have now expressed both sides of Eq. (43) through the coordinate time t' a . Since Eq. 
(44) is valid at any time, t' a is a mere "dummy" variable; it may be denoted by t again or 
even be dropped as a function argument. This yields the following expression for the action 
of the generator Gi on the spatial coordinates x valid for any boost vector Vi 

J2 V i \ h X ajXai ~ SijtJ = 5 v X ai = J2vj{x ai , Gj}* + O ( J , (45) 

with n being again the order of approximation. Inserting the equations of motion (28) into 
Eq. (45), we deduce the world line condition 



{x a i,Gj} — ^ 2 X a j{x a i,H} 5ijt + O ^^(n+l) ^ ' 



(46) 



This condition joined to the Poincare algebra, to the expression of the action of Pi, Ji on 
the phase space variables, and to the physical requirement that H is the Hamiltonian of 
the system will be sufficient to determine Gi uniquely (up to arbitrary generalized canonical 
transformations) . 

In the following, it will be useful to perform computations with the quantity K; L = Gi + P{t 
instead of Gi. Because of Eqs. (35a, 35b, 35d, 35i'), Ki fulfills the same Dirac-bracket 
relations as Gi, but is not explicitly time dependent; it is not an integral of motion either. 
Indeed, we conclude from Eq. (35g) combined with the conservation law ^ = that 
^ = -Pi. Applying the definition of K t we find ^ = and {K^H}* = P^ The 
world-line condition for Kj reads, 



1 / 1 

{X a i, Kj} = -^X a j{x ai , H} + O ( c2 ( n+1 ) 



(47) 



It can be derived from Eq. (46) and from the first of the two equations (38a), which, due to 
the first class property of Pi and relation (27), also holds as Dirac-bracket relation on the 
constraint surface as mentioned above. 

V. APPLICATION TO 3PN DYNAMICS 

The conservative part of 3pN equations of motion for compact binaries in harmonic 
coordinates allows for a Lagrangian of the form [21] 

2 

2 



1 2 1 

L(X,X (1) ,X (2) ) =-J2 m a(^a ) ) 2 + V (x) + -V 1 (x,X^) 



a=l 
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+ ±V 2 (x,xW,xW) + ±V 3 (x, x«,x (2) ) + (48) 

with Vo(x) = G mi ™ 2 . Beyond the fact it is restricted to a two-body system, it does not 
show any dependence on x® in the term V3 by contrast to the general Lagrangian (6). 
The investigation of the dynamics can be performed in two different ways: First, since the 
original formalism as displayed in sections II and III does not require the invertibility of 
the Hessian, it may be applied on an artificially extended configuration space putting up 
however with lengthier computations. Second, we may consider V%{x, x^\ x^) as a higher 
order correction to V 2 (x, x^), gathering both functions into a new quantity, 

V 2 (x : x^\ x^) = V 2 (x, x W, x W) + \v 3 (x,xW,x<V). (49) 

cr 

In this approach we write the 3pN Lagrangian (substituting e = ^) as 
1 2 

^3 P n =2 E m «(4 1) ) 2 + V (x) + eV^x, x«) + e 2 V 2 (x, x«, x^) + O (e 4 ) , (50) 

a=l 

without affecting the original accuracy. The formalism described in sections II and III has 
to be appropriately modified. 

It has been checked explicitly that both methods lead to equivalent constraints, to iden- 
tical expressions for the elementary Dirac-brackets and to identical conserved quantities on 
the constraint surface. However, since the second method is more adapted to the problem, 
it demands less computational effort. Furthermore it can be employed for the computations 
of 2pC-dynamics when the terms of third order are neglected. We shall therefore use it for 
the subsequent calculations. 

The equations of motion for the Lagrangian (50) read 

2 

-m a xW +J2e s A as (x,...,x^) = 0(e 4 ), (51) 

s=0 

dV dV L _d_dV L dV^ _ d_dV^ d 2 dV 2 

a0 ~dx~ a ' al ~^'^^' Aa2 -dx~ a ~Jt^m + dT^- 

The computation of the constraints for a system described by the above Lagrangian is 
performed in appendix B. The main difference to the formalism exposed in Sec. II is that 
the equations of motion have to be employed twice at the end to guarantee time stability. 

Finally we obtain a minimal set of time stable Lagrangian constraints (B7) 

2 

m a xt ] - 5> S £<w(x,x«) = 0{e A ), (52) 
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for a = ai and r = 2, 3. The Ostrogradski momenta derived from (29) are given by 

n 0a = m a x^ + e<5> 0a {x, x (3) ) + 0{e 4 ), (53a) 
n la = e 2 <Z> la (x, X W,xW) + 0(e 4 ), (53b) 

where the $o Q , ®ia are now constructed with V\ and V^- 

dV 1 dV 2 d dV 2 , x 

^oa = — 777 + £ — 777 — £^ 757 (54a) 

<9a4 J 8x a ' dt Q Xa ' 

dVo 

*'» = SF- (54b) 

The transformation of the Lagrangian constraints into constraints on the phase space coor- 
dinates yields the Hamiltonian constraints 

Xi a = n lQ - e 2 <£ la (x, x&) = 0(e A ), (55a) 

^ a = x« - — [n 0a - e$ 0a (x,xW)} = 0(e 4 ). (55b) 
m a 

The $j a (x, x^) J s are derived from the $ JQ , (rr, a^ 3- -^) 's by eliminating higher order deriva- 
tives by means of Eqs. (52). The Poisson-bracket matrix D of the constraints (55a, 55b) is 
regular, as shown in appendix C by an explicit computation of D~ l . It is thus possible to 
endow phase space with the Dirac brackets (25). Their explicit expression can be found at 
the end of appendix C. 

We are now in position to give the representation of the generators of the Poincare group 
as generators of infinitesimal generalized canonical transformations with respect to the Dirac 
bracket for the conservative 3pN binary dynamics in harmonic coordinates. The phase space 
variables X\ and x 2 are the positions of the two point masses. The generators act as usual 
on them and on their conjugate momenta, in the way specified in Sec. IV. 

The reduced Hamiltonian, generator of infinitesimal time translation, can be computed 
from 

2 3 2 3 

#3pN = -^3pN + 110 aiX ^i + YYl Ul aixi S ' ( 56 ) 

a=l i=l a=l i=l 

by using the constraints to eliminate x^\ x^ and Hi in favor of the coordinates of the 
constraint surface. The latter are chosen to be x and Il because of their "Hamiltonian" 
character (a possible alternative is x and x^). In this grid, the reduced 3pN Hamiltonian 
explicitly reads 

#r,3pN = °#N + + + — 3 H N , (57) 
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We have again adapted the notation of Ref. [21] introduced after Eq. (5). Furthermore, we 
have posed n i2 = Xl ~ X2 . The term "+1 < — > 2" represents the expression that precedes but 
with interchanged particle indices, including the contributions that are already symmetric 
in the particle indices, particularly ni 2 must be changed into n 2 i = — riu there. Note that 
it is possible to remove the logarithm terms through a coordinate transformation preserving 
the harmonicity conditions outside the bodies [10]. 

The generators of infinitesimal spatial translations and rotations are given on the whole 
phase space by Eqs. (40, 41) as 

2 

Pi = ^^n 0a i = n ii + n 02 j, (58) 



a=l 
12 3 



T - V" V" r -r (m) TT 

"i / j / / ^ijk-^aj 11 mak 

m=0 a=l j,k=l 
2 3 2 3 

= ^ ^ ^ijkXaj^-Oak + ^ ^ tijkXaj^-lak- (59) 
o=l j,k=l a=l j,k=l 
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For the dynamics, only their restrictions to the constraint surface are relevant as integrals of 
motion and generators of symmetry transformations. They can be computed by eliminating 
and 111 from Eqs. (58, 59) with the help of the constraints (55a, 55b). The result is 



Pi\r — rioii + n 2i, 

2 3 

^Ir = €i i k 

a=l j,k=l 
1 
TO, 



(60) 
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aj 



■e 2 — <f> laj {x, x {1) (x, n )) j U 0ak 
m a J 



(61) 



-e 3 — ® 0aj (x,x (1 \x, Tl ))® lak (x,x {1) (x,Tlo)) 
The angular momentum Jj| r of the post-Newtonian dynamics can be written explicitly as 

(62) 
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where a sum over the indices j and k = 1,2,3 must be understood. It has been checked 
that the above momenta satisfy the relations (35a-35e). In particular, H is invariant under 
spatial translations and rotations and the components of the generalized total linear and 
angular momentum, and Jj, are integrals of motion. 

The generator (62) takes a simple familiar form when written in coordinates that are 
canonically conjugate with respect to the Dirac bracket. We have obtained these coordinates, 
expressed in terms of x, n by guess work. They read 

P a = n 0a - e 3 Ho) d% f' U ° ] + O ( g 4 ) , (63a) 

m 7 ox a v ' 



Q a = x a - e 2 —$ la (x, Ho) +e 3 J2 — <M*. n )^p^ + O (e 4 ) , (63b) 
m a m^ OL\.Q a 

where the sum holds over all pairs of indices 7 = bj, b = 1,2, j = 1,2,3 and where 
<3> S7 (a;, n ), s = 0, 1, is a short notation for & sy (x, x^(x, n )). This result can be verified by 
using the explicit expressions (C12-C14) for the elementary Dirac brackets. 

It is important to stress that there does not exist any harmonic coordinate system in 
which Q and more generally any other canonical coordinates represent the particle posi- 
tions. If it were possible, we could go over to a frame where the particle positions are 
canonical coordinates while maintaining manifest-Poincare invariance. In terms of these 
new coordinates, the Hamiltonian as well as the Lagrangian would be ordinary although 
describing manifest Poincare invariant dynamics, in contradiction to the no-interaction the- 
orem [24]. A consequence is that we are not able to obtain canonically conjugate particle 
coordinates by means of Poincare transformations. This can be seen from Eqs. (C15, C16), 
which imply that {x a ,xp}* cannot vanish everywhere for the 3pN and 2pC dynamics. The 
form of these two relations is maintained by Poincare transformations that are also gener- 

25 



alized canonical transformations. However, in the case of vanishing coupling constant, we 
have Q — > x, P — > n . In other words, for the non-interacting systems, manifest Poincare 
invariance becomes compatible with the choice of spatial coordinates as canonical coordi- 
nates. 

The interest of using Q and P rather than x and n for explicit calculations is that in 
these coordinates the Dirac bracket takes the simple standard form of a Poisson bracket. 
For simplicity, we shall note any function F(x, Ho) on the constraint surface expressed by 
means of the canonically conjugate variables Q, P as F = F(Q, P) = F(x(Q, P), Ho(Q, P))- 
The Dirac bracket at a point x(Q, P),U (Q, P) of the constraint surface then reduces to 



{F,GY 



r ~ ~. dF dG dG OF 

= {F, G} Q , P = ^ QQ-Qp -QQ-Qp- ( 64 ) 
x(Q,P),Tl (Q,P) a a Va ur a y V« ur ot 



In canonically conjugate coordinates, Pj and Jj have their usual expression on the con- 
straint surface for both dynamics under consideration, namely 

2 

Pi = ]TP« + 0(^ +1 ), (65) 
0=1 
2 3 

Ji = Yl Yl e ^Qa 3 Pak + O {e n+1 ) , (66) 

a=l j,k=l 

with n = 3 at 3pN (n — 2 at 2pC). Indeed, spatial translations and rotations leave the 
constraint surface invariant so that the restrictions of Pj and Jj to V are the well known 
generators of infinitesimal spatial translations and rotations on Y (cf. Sec. V) taking above 
shape in terms of canonically conjugate coordinates. This explains incidentally the absence 
of a first post-Newtonian contribution in the expressions (62) and (83) for Jj| r as a function 
of x, n . If the term containing IT does not contribute to the first order, this is because the 
Ill's are second order quantities; for the rest, the x's and the n 's are canonically conjugate 
modulo 1/c 4 corrections. 

To determine the generator of Lorentz boosts, it is useful to work with the canonical 
coordinates Q, P on the constraint surface. In addition, it is more convenient to derive 
Ki = Ki(Q,P) = Gi(Q,P) + Pi(Q,P)t rather than G;(Q,P) itself. We compute k t with 
the help of the method of undetermined coefficients [17], thereafter fixing the uniqueness of 
the solution. 

As already mentioned, Ki fulfills the same Dirac bracket relations (35f-35h, 35i') as Gj, 
but it is not explicitly time dependent and it is not an integral of motion. The behavior of 
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Ki under spatial rotations is governed by Eq. (35h) with the angular momentum Jj given 
by (66). We conclude that Ki has the general structure 

2 

£ = £ M a (Q, P)Q m + N a (Q, P)P ai + O (e n+l ) , (67) 

a=l 

where M a (Q,P), N a (Q,P) are two post-Newtonian scalar functions. The form of the dif- 
ferential equations for Ki resulting from the Poincare algebra and the world-line condition 
suggests for M a (Q,P) and N a (Q,P) the ansatz 

Cn ,...,n 5 R no Pi ni Pt 2 (PiP2Y ls (N 12 P 1 ) n ^N 12 P 2 y^ + d 1>Sl R s ^ ]n(R/n) + d 2yS2 R s > \n(R/r 2 ), 

(68) 

with R — \Qi — Q 2 \ and Ni 2 = Ql ~^ 2 ; the logarithm terms are only expected to appear in 
the function M a and the powers n ,si,s 2 are presumed to be integers, while ni,...,n 5 are 
natural numbers. The admissible eight-tuples of s±, s 2 , n , n 5 are restricted by demanding 
the correct physical dimension for M a and iV a . 

Insertion of the preceding ansatz into the partial differential equations (35f-35h, 35i') 
yields the linear system of equations to be solved. The coefficients seem to be over- 
determined, but the equations are not all independent so that there actually exists a so- 
lution. It is most easily derived by means of a Computer-algebra program 1 . If we only 
impose that the Poincare algebra should be satisfied, some coefficients remain undetermined 
even at the lpN order. They are set by requiring the world-line condition on Ki, reflecting 
the manifest Lorentz invariance of the system. As a result, the expression for Ki obtained 
through this procedure, is automatically consistent with both the Poincare-algebra and the 
world-line condition. 

Now, as discussed before, the canonically conjugate coordinates simplifying our explicit 
calculation are not harmonic. Re-expressing Ki by the coordinates x and n of the constraint 
surface we finally arrive at Ki in a harmonic-coordinate frame. The result is displayed up 
to third post-Newtonian order as 

K^p^i = K m + — K m + — K Ni + — K Ni , (69) 

°K Ni =m l xu + m 2 x 2i , 
1 We have used the software Maple 8 ©Waterloo Maple Inc. 



27 



n 2 



n r 2 



2mi 



Oj Ll()2 Gil),!!),. 

Xli + x 2i {x u + x 2i ), 



2m 9 



2r 



n i 4 g 



8mi 3 Li ' 8m 1 m 2 
- 14m 2 2 (ni 2 n i)n ij 

G ?2 mim2 



7m 2 2 noi 2 ni2i - m 2 2 (ni 2 II 01 ) 2 n 12i - 14mim 2 (n 12 n 02 ) n ii 



G 

+ — 

4r 



+ 7 (n iIIo 2 ) xu 
G 



+ 



4 r i 



n i 2 xii + (ni 2 n i) (ni 2 n 02 ) x u 

+ 1< — >2, 



mi 

— 5mia;H + 7m 2 xi 



n 6 



42m 2 3 „ , 



6m 



mi 



n i V 2i — (ni 2 n i) n i 2 n 



mi 



16mi 5 48mi 2 m 2 2 

- 6m 2 2 (ni 2 Il 02 ) n i 2 II ij - 30mim 2 (ni 2 II 02 ) 2 II i 2 ni 2 ; 

- 24m 2 2 (ni 2 n 02 ) (ni 2 n i) n i 2 ni 2i - (ni 2 n i) 2 n i 2 ni 2i 

mi 



Oli 



+ 18mim 2 (ni 2 n 02 ) 2 (n^IIoi) n ii - 48mim 2 II 02 2 (n^IIoi) n 
+ 12mim 2 (ni 2 n 02 ) (n ill 02 ) U 01i + 20mi 2 (ni 2 II 02 ) 3 II 



Oli 



L 01i 



+ 



20m 2 3 
mi 



(ni 2 noi) 3 n ii + 18m 2 2 (n 12 n 02 ) (ni 2 n i) 2 n ii 



- 48m! 2 n 02 2 (ni 2 n 02 ) n i, + i2m 2 2 (n 12 n i) (n in 02 ) n 

(ni 2 n i) 4 ni 2i + 6m 2 2 (ni 2 n 02 ) (n 12 U 01 f n 12i 



Oli 



+ 



3m 2 3 
mi 



+ 18m 2 2 ni 2i (n 12 n i) 2 (n iIIo 2 ) 



+ 



G 

16 r 



'21m 2 



14 



3 IIoi 4 a;H + ( n i2n 02 ) 2 n 01 2 a;ii + — - (n in 02 ) n i 2 a;ij 



mf m\m 2 

'2 



mi" 



2 

-n 02 2 n i 2 xii (ni 2 n i) (ni 2 n 02 ) iioi 2 ^ (n in 02 ) 2 a;H 

m\m 2 mf m\m 2 

12 2 

(ni 2 n i) (n 12 n 02 ) (n iiio 2 ) x u h — -n 02 2 (n in 02 )x H 

m\m 2 m 2 l 

2 q ^22 
r n 02 (ni 2 n 02 ) (ni 2 n i)a;ii (ni 2 n i) (ni 2 n 02 ) 



m 2 " 
llmi n 4 r 

m 2 3 



m\m 2 



+ 



2 68m 2 - ., > l//; 2 

- 299m 2 n i n 12i H n i n i2i + 



34m- 2 



mi 



mi 



(ni 2 n i) n 



Oli 



48 r 

+ 428m 2 (ni 2 II i) n ii + 386m 2 (n 12 II 02 ) n ij + 76mi (ni 2 n 02 ) II 

+ 48mi (ni 2 II 02 ) (ni 2 n i) n 12i - 100m2 (ni 2 n i) 2 n 12i + 13m 2 (ni 2 n i) 2 n 12i 

mi 



Oli 



28 



+ 



+ 816toi (n in 02 )ni2 

1 G 2 



3360 r 2 
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The 3pN accurate generators H 3pN ,P 3pNi , J3 P Ni,G 3pNi = K 3pNi - P 3pNi t we have derived 
in harmonic coordinates, constitute an approximate representation of the Poincare Algebra 
given by the generators of infinitesimal canonical transformations with respect to the Dirac 
bracket, and thus, they generate the Poincare transformation group on the phase space 
coordinates. The quantities are integrals of motion that may be interpreted physically as 
the 3pN conserved total energy, the generalized total linear momentum, the generalized total 
angular momentum and the center of mass constant. Their expressions have been checked 
up to first order by comparing them with Ref. [43]. If we adapt the "Lagrangian-like" 
coordinates x, x^ on the constraint surface, the above generators of infinitesimal canonical 
transformations representing conserved quantities reproduce the Noetherian constants of 
motion associated with the Poincare symmetry of the dynamics that are derived in [21] by 
means of the Lagrangian formalism. Note that the definitions of Ki and Gi in this article 
are reverse. 

It remains to show the uniqueness of the generators Pi, Ji, Gi of the Poincare group up 
to generalized canonical transformations provided that if is a post-Newtonian Hamiltonian. 
Let us indicate that the following proof is also valid for the post-Coulombian dynamics. 

We start with the case of the momenta Pi and Jj. Their uniqueness as a function of 
the phase space coordinates has been established on the whole phase space in Sec. (IV). It 
entails the uniqueness of their restriction to the constraint surface. 

We next consider the center of mass constant Gj. Its uniqueness is equivalent to that 
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of Ki = Ki(Q,P) = Gi(Q,P) + Pi(Q,P)t, since the coordinate transformation is a dif- 
feomorphism and since p is unique. Again, Q, P denotes the set of canonically conjugate 
coordinates on the constraint surface while p represents the total linear momentum of the 
system. Let us assume that there exist two solutions Ki and K[ = Ki + fi(Q, P) that both 
approximately fulfill the Poincare algebra and the world-line condition within the common 
coordinate frame. The relations (46, 35g, 35h) written in terms of the coordinates Q, P as 
given by Eqs. (63a, 63b) read (cf. also (64)): 



Q,p 



+ O (e 4 ) = ex aj {x ai , H} QtP + O (e 4 ) , (70) 
{K j ,H} QtP =P j + 0{e'), (71) 

3 

{Ji, K 3 } Q , P = tijkKk + O (e 4 ) . (72) 
fc=i 

Note that, because x a , Hop and Q a , Pp agree at zeroth order, the Newtonian (or Coulombian) 
contribution to the Hamiltonian in canonically conjugate coordinates has the form 

°H=Y,^+ u (Qy ( 73 ) 

a=l a 

The uniqueness of Ki is then proved order by order using Eqs. (70-73). First, we may insert 
successively Kj and K'j into (70) since they are both assumed to fulfill the latter relation. 
Taking the zeroth order of the difference, we find that 

d%(Q,P) 



dp 



= 0, (74) 



at 



hence °fj(Q,P) = °fj(Q). Next, we go over to Eq. (71), where we insert again Kj and K'- 
successively before subtracting the ensuant equalities. By virtue of Eq. (73) we obtain at 
zeroth order: 



3 



a=l k=l 

Because Q a ,P/3 constitute a set of independent coordinates, and Eq. (75) holds for all P a k, 
we are led to 



d%(Q) 
dQ 



ah: 



%(Q) = const. (76) 
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From Eq. (72) taken at zeroth order, we conclude by applying the same procedure as ex- 
plained above, that 



= {%%} q ,p = {%°K' 1 -%} q ,p 

3 3 
k=l k=l 

=► % = =► K: = °K t . (77) 



The proof of uniqueness at the first order is similar. The world-line condition Eq. (70), 
truncated at this level of approximation yields 1 fj{Q,P) = 1 fj{Q). Eq. (71) reduces to 

{%,°H} Q , P + {%, 1 H} Q , P = 0, (78) 

and we have the same equation for K'y On the other hand, we know from Eq. (77) that K[ 
and Ki differ at most from the first order on. Therefore taking the difference between the 
equation (78) for K'j and the same equation for Kj at first order leads to 

CfjiQ), °H}q,p = 

2 3 



Pgkd l fj(Q) = Q 

^-i ^ m a dQ ak 

a=l k=l ^ 



d'fAQ) 



0, (79) 



9Qak 

so that 1 fj{Q) = const. Now, from the first order truncated version of Eq. (72), it follows 



= {% 1 f 3 } Q , P = {% 1 K>- %} Q:P 

3 3 



k=l k=l 

=S> % = x k[ = %. (80) 

The uniqueness of second, third or even higher orders follows analogously. 

VI. APPLICATION TO 2PC DYNAMICS 

The Lagrangian describing the binary dynamics of the Feynman- Wheeler theory in 
Lorentz gauge up to second post-Coulombian order takes the form 



L(x,x (1) ,x (2) ) =^m a Or«) 2 + \/ (x) + ^Vfax™) + \v 2 (x, x (2) ), (81) 



a=l 
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with Vo(x) = —^f 1 . We may thus adopt directly all general results derived for the 3pN 
dynamics by neglecting the third order contribution. The elementary Dirac brackets of the 
2pC dynamics, for instance, may be inferred from appendix C. The computation of the 
generators of the Poincare group corresponding to the conserved quantities of the Feynman- 
Wheeler 2pC binary dynamics can be performed in a way similar to the post-Newtonian 
case. The results read as follows. 

The reduced Hamiltonian is given by 



I i 

#2 P c = H c + -g H c + -r H c , 



(82) 
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This result confirms the one derived in Ref. [34] as corrected in Ref. [35]. 

The generators of spatial translations and rotations Pi and J« can be derived from Eqs. 
(60) and (61). The explicit expressions for the restrictions on the constraint surface in terms 
of x and II read 



— n ii + n 02 i, 

^ 2 pCi = °Jd + — 2 Jd, 



(83) 
(84) 



ijk-^lj no 2 fc, 
eie 2 



2 T —c 

u Ci — t ijk 



! ""'4m2mi 2 
e eie 2 
yfc 8m 1 2 m 2 2 r 



rioi^i 2 ) n ijn 02 fe + 



eie 2 



8mi 2 m 2 2 r 

2 



m 2 n i + m 2 (n ini 2 ) 



miil 02 2 - m 1 (n 02 ni 2 ) z 

^ljrTlo 2 A; 



xi 7 -n, 



ei 2 e 2 2 
4mim 2 r 2 



XiTI 



ljJ-J-oifc — ^ijn 02 A; 



+ 1 < — ► 2, 
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where a sum over the indices j and k — 1,2,3 must be understood. 

The generator Gi is computed as in the post-Newtonian case. For Ki = d + Pit in terms 
of x and n , we have 



°K Ci =m l x li + m 2 x 2 i, 

x v n i 2 n 2 2 . eie 2 

K Ci =7, X U + 7, X 2i + T. [ x li + x 2i\, 

2m i 2m 2 2r 



(85) 



2 K Cl = 



^l^noi 4 



eie 2 



+ 



ami 
2m 2 
mi 



8mim 2 

(ni 2 n i) n ij 

n\2i + 1 < — 



(ni 2 n 01 ) niM n i 2 n 12i + 2 (ni 2 n 02 ) n ii 

mi mi 



eie 2 



4mim 2 r 



^n iii 02 ) xii + (ni 2 n i) (ni 2 n 02 ) x u 



2. 



4m 2 r 

The above generators H, P iy Jj and Gj of infinitesimal generalized canonical transformations 
provide a representation of the Poincare algebra on the constraint surface of 2pC Feynman- 
Wheeler binary dynamics with respect to the Dirac bracket. The uniqueness of Pj, Jj and 
Gi up to a generalized canonical transformation has been established in Sec. IV. 



VII. SUMMARY AND DISCUSSION 



In the present article, we applied an appropriate canonical formalism to describe the 
third post-Newtonian dynamics of point-mass binaries in harmonic coordinates. We treated 
the second order post-Coulombian dynamics of Feynman- Wheeler theory in Lorentz gauge 
analogously. In contrast to earlier works, we did not leave the coordinate conditions by 
performing a higher order contact transformation [30] but instead we generalized a method 
developed in Ref. [34] and constructed the dynamics directly in harmonic coordinates resp. 
Lorentz gauge within the framework of canonical formalism both singular and of higher order 
in the time derivatives. The canonical formulation opens the way to advanced investigations 
about the geometrical and physical interpretation of the motion. It is highly desirable 
for deriving the generators of infinitesimal generalized canonical symmetry transformations 
that provide the integrals of motion. We computed for the first time the generators of the 
Poincare transformation group or, equivalently, the conserved quantities corresponding to 
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the manifest Poincare invariance, for third post-Newtonian conservative binary dynamics 
in harmonic coordinates as well as for Feynman- Wheeler second post-Coulombian binary 
dynamics in Lorentz gauge. An appropriate choice of coordinates of the constraint surface 
reveals that the 3pN conserved quantities we have obtained agree with those derived via the 
generalized Noether theorem in [21]. After being reduced to the center of mass frame [40], 
they can be used for the derivation of an analytic parametric solution to the third post- 
Newtonian equations of motion in harmonic gauge for compact binaries in eccentric orbits 
[41], which is in turn of high practical relevance for the construction of gravitational wave 
search templates and comparison with numerical simulations. A useful generalization of the 
present investigations is suggested to be the application of our method to post-Newtonian 
binary dynamics in harmonic coordinates including spin. 
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Appendix A: LAGRANGIAN CONSTRAINTS [34] 

In this appendix we display for completeness the iterative method developed in Ref. [34] 
to derive a minimal set of time stable constraints for a dynamics described by a Lagrangian 
with the structure shown in Eq. (6). 

(i) The first stage of the method consists in transforming the primary constraints into 
a more restrictive set implied by the primary constraints through time stability, without 
making use of the time evolution equations. In this new set, the higher order time derivatives 
are found separately in the various equations. 

The primary constraints (10) imply 

e^CO = £ nJ- Aa0 (x) + (e n+1 ). (Al) 
m a 

Thus, Eqs. (10) multiplied by e n ~ x are fulfilled by the motion, as well as the relations 



34 



following from (Al) through repeated time differentiation 



£ n x (2 + r) = £ n^ Bat2+rfi ( x , x W) + 0( £ n+1 ), 



rn 



(A2) 



where r = 0, ...,2n — 3; The quantities B at 2+ r ,o(x, x^) are derived from the A q0 's by differ- 
entiating r times and eliminating accelerations by virtue of Eq. (Al) whenever they occur; 
the relation (A2) reduces to (Al) when r = 0. Beware that, because of the occurrence of the 
time derivative Xa n \ Eq. (A2) with r = 2n — 2 is not a constraint, but rather the condition 
of time stability for the actual constraints given by the system (A2) for r = 0, 2n — 3. It 
becomes a consequence of the equations of motion as soon as the dynamics is restricted to 
the constraint surface. 

The next set of constraints implied by the primary ones can be found by multiplying Eqs. 
(10) by e n ~ 2 : 

1 



[A a0 (x)+eA al (x,x^,x^)]+O(e n 



(A3) 



The acceleration dependence of A al (x,x^\x^) may be eliminated by means of Eqs. (A2) 



£ n-l x (2) 



.n-1 1 



[A a0 (x)+eB at2<1 (x,x^)] +0(e n+1 ). 



(A4) 



Differentiating r = 0, ...,2n — 3 times with respect to time and replacing all occurring 
accelerations by means of Eqs. (A4) results in 



e »-i x (,2+r) = e n-iJ_ [B a , 2+rfi (x,x^) + eB a , 2+rA (x,x^)] +0(e n+1 ). 

TTh, 



(A5) 



Carrying on with this procedure we arrive at the set of constraints 

"n-1 



ex% +r) = e- 



^2e s B a ^ s {x,x (1) ) 



s=0 



(A6) 



where r — 0, 2n — 3. 

(ii) In the second part of the method proposed by Jaen, Llosa and Molina, the above 
system is completed to one that is stable under the time evolution ruled by the Euler- 
Lagrange equations. The time stability condition for Eqs. (A6) which also covers the time 
stability of all previous sets of constraints, equals (A6) with r = 2n — 2: 

"n-1 



sx^ = e 



a,2n,s{X, 



X 



s=0 



+ 0(e 



n+l\ 



(A7) 
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Demanding time stability of the constraints for the motion on the constraint surface is 
equivalent to removing higher order time derivatives in the Euler-Lagrange equations with 
the help of Eqs. (A6) and eliminating Xa^ by means of Eq. (A7). We obtain the new set of 
constraints 



x 



(2+r) 



1 



(1)^ 



+ 0(e n+l ), 



(A8) 



with r = 0, 2n — 3. From the way this set of constraints has been obtained, we see that 
it must hold even if some of the Xa^ do not occur in certain linear combinations of the 



Euler-Lagrange equations, which may be the case for non-invertible 
For the set of constraints (A8), the time stability condition reads 



d 2 Vn 
OX ai OX bj 



T (2n) _ 



+ 0(e n+l ) 



(A9) 



This condition is satisfied in the sense that removing the higher derivatives in the Euler- 
Lagrange equations with the help of Eqs. (A8) and substituting xi 2n ' by Eq. (A9) yields a 
constraint (Eq. (A8) with r = 0) fulfilled on the constraint surface. It is obvious, that the 
set (A8) is more restrictive than the original primary constraints. Nonetheless, we see from 
the derivation that it follows from them through the time stability condition. As a result, 
with (A8), we have found a minimal set of constraints stable in time for the dynamics under 
investigation. 



Appendix B: 3PN CONSTRAINTS 

In this appendix we derive the constraints for the two-body system at the third post 
Newtonian order by modifying appropriately the method displayed in appendix A. 

Within our approximation scheme, any vector being a multiple of e 2 is a null- vector of the 
Hessian of the 3pN Lagrangian (50). Therefore, using Eq. (51) as initial primary constraints, 
we have 

e 2 (m a x& - ^e s A as (x,...,x^ = 0(e 4 ). (Bl) 

Primary constraints due to the possible existence of further null vectors are not precluded. 

(i) We now transform the system (Bl) into a more restrictive set obtained by implying 
time stability, without employing the explicit equations of motion yet. 



36 



The acceleration dependences in orders higher than e 2 can be removed from Eq. (Bl) by 
means of the constraints obtained through multiplying the primary constraints with e. We 
are led to 

e 2 (m a x® -J2e s B a , 2tS (x,x^)j = 0{e A ). (B2) 

After differentiating and eliminating the occurring accelerations by making use of Eqs. (B2), 
we find 

e 2 (max® - ^e-B a , 3 ,.(x,x«) j = 0(e 4 ). (B3) 

Requiring that the above constraints are satisfied for all times yields the time stability 
condition 

e 2 (m a x® -J2e s B aAs (x,x^)j = 0(e 4 ). (B4) 

(ii) Additional constraints emerge from the latter time stability condition as well as from 
the equations of motion. They are derived by multiplying the Euler-Lagrange equations by 
e and thereafter eliminating all second and higher order time derivatives and found beyond 
the leading term with the help of Eqs. (B2-B4). A last time differentiation and elimination 
of the newly appeared accelerations leads to the constraints 

e (m^ -J^e s B a , r , s (x,xW)j = 0(5 4 ), r = 2,3, (B5) 

which are not yet time stable. The time stability condition 

e (m a x^ -J2e s B aAs (x,x^)j = 0(e 4 ) (B6) 

inserted together with (B5) into the equations of motion yields the more restrictive set of 
constraints 

2 

™ a x^ -^e B B a<r>8 (x,xW) = 0(e A ), r = 2, 3. (B7) 

s=0 

These constraints are stable, for inserting all constraints plus the time stability condition 
into the equations of motion yields a constraint (namely (B7) with r = 2) that is actually 
satisfied on the constraint surface. In contrast to appendix A, the time stability condition, 
and thus the equations of motion had to be employed twice to determine the constraints of 
3pN dynamics. Finally, we note that it is now possible to argue as in Sec. II that further 
null vectors of the Hessian do not lead to further constraints. 
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Appendix C: 3PN ELEMENTARY DIRAC BRACKETS 



The issue of this appendix is to present the explicit computation of the Dirac brackets 
for the dynamics of interest in this article. The Dirac bracket is bilinear and acts as a 
derivation on each argument. To know its action on any two functions of the phase space, 
it is therefore sufficient to know all elementary Dirac brackets, i.e. the Dirac brackets of the 
constraint surface coordinates. In order to compute them, we have first to determine the 
Poisson-bracket matrix D of the constraints up to the needed third order in e, appropriately 
adapting the general method given by [34] to the 3pN case. All expressions reduced to 
second order in e can be directly used for the post-Coulombian dynamics case. 

We start by splitting D into 6 x 6-submatrices as 

I s t\ 

D=\ , (CI) 

\-T T U 



So* = {Xia, Xip} = - e 2 ^% + O (e 4 ) , (C2) 

OXa OXp 



= (Xia, ^} = -8* - e—^ + e —— + 0(6), (C3) 

U a , p = {coia, = e^— - ^) + O (e 4 ) . (C4) 

m a mf3 \ dx a 0Xf3 J 

This matrix is clearly invertible. To obtain iteratively the inverse matrix of D, we write the 
expansion in powers of e as 

n 

d = j2 £SSD +°( £n+1 )- ( C5 ) 

It can be easily verified, that D~ l up to n-th order is given by 

D~ x = ^l + ^e«wj D- 1 + O(e n+1 ), (C6) 
where S N is computed iteratively from 

s-1 

S N=- ^D- 1 S D-J2 ® D ~ l S ~ rD rjv - ( C7 ) 

r=l 

In our problem, we need D~ x up to the third order. Its explicit expression reads 
D' 1 = (l + J2 £SSN ) D- 1 + O(e 4 ) 



s=l 



38 



= D -l _ £ D -1 1 D Q D -1 +e 2 ( _0^-l 2 D Q D -1 + D -1 1 D D -1 1 D Op-l) 
+ £ 3 (_ D -1 S D D -1 + D -1 2 D D -1 1 D D -1 + D -1 1 D D -1 2 D Op-l 

_0 D -1 1 D D -1 1 D D -1 1 D D -lj +O ( £ 4y ( C g) 



At last, we split D 1 into submatrices, 

D- l = 




(C9) 



and account for 



{U 0a , uip} 
{n 0Q , Xip] 



0(e°) , 

0(6 4 ) , 

0(e), 



From the expression (33) for the Dirac bracket we conclude directly that, to determine the 
elementary Dirac brackets up to the desired third order, we need Y up to order 1, Z up to 
order 2 and that we do not need X at all. The components of the relevant submatrices are 







1 


dxi 1 ^ 


—5 a/3 


— e — 

m a 




=e 2 




dx^ 














\dxf 


dx? 



(CIO) 



E 



(9$ 



1 7 



7 

1 (9$ 



dx 



(i) 



ggla 



1 9$ 



0-> 



™ 7 ftr^ 



07 



™ 7 cfei 1} 



0(e 4 ). 



(Cll) 



Hence the elementary Dirac brackets read 
1 _ 



{x a ,x p }* -- 

m a mp 

{x a ,Up}* =5 al3 + e 



J a,f3 



(C12) 



2 1 d$> la 



m a dxp 

' J m . m 



«9$ la \<9$ 07 <9$ 07 <9$ 



dx 



(i) 



o 7 



/ dx P dx ( a ] dx p 



+ 0{e 4 ) , 



^ m 7 \ ox a 



Q 7 U^ly 



dx 



a 



<9$o 7 9$i 7 



0{e% 



(C13) 
(C14) 
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Finally, we explicitly display the relation (C12) up to second order for both 3pN and 2pC 
dynamics. 



G 

2pN: {x ai ,x bj }* =e 2 — 



7 



n 



ai 



n 



hi 



m n 



m b 



n 



Obk 



Wabi 



+ (7dij + n abi n abj ) 2_^n abk I 1 



(C15) 



2pC: {x ai ,x bj }* = - e 



2 ^a^b 

m a m b 



1 n a6 j - 

m a m b ■ 

3 



n 



TO 6 



m a 



Tt'ab i 



+ (Sij + n abi n abj ) ( ^Tn abk | + Ih 
\k=i \ ma m 



(C16) 
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